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Abstract 



Rigorous justification of the Hubbard-Stratonovich transformation for the Pruisken-Schafer 
type of parameterisations of real hyperbolic 0(m, n)— invariant domains remains a challeng- 
ing problem. We show that a naive choice of the volume clement invalidates the trans- 
formation, and put forward a conjecture about the correct form which ensures the desired 
structure. The conjecture is supported by complete analytic solution of the problem for 
groups 0(1, 1) and 0(2, 1), and by a method combining analytical calculations with a simple 
numerical evaluation of a two-dimensional integral in the case of the group 0(2, 2). 

1 Introduction and formulation of the conjecture 

For more than two decades, the nonlinear er-model methodology has been widely applied to 
studies of single electron motions in disordered and chaotic mesoscopic systems [TJ [5]. The 
method was pioneered by Wegner [3j and further developed by Wegner and Schafer [3] , and 
Pruisken and Schafer [5] in the framework of the replica method used to reduce one-particle 
Hamiltonians with microscopic disorder to a nonlinear er-model. In the early eighties, Efe- 
tov [6] introduced the supersymmetric variant of the method which avoided the problematic 
replica trick and directly led to the supermatrix version of the nonlinear er-model. Since then 
this latter nonlinear er-model has been also successfully applied to a variety of problems in 
the framework of random matrix approach to chaotic scattering [7] [8], Quantum Chromo- 
dynamics [5], as well as a few other fields of physics. 

A standard derivation of the nonlinear er-models requires to use at some point the so- 
called Hubbard-Stratonovich transformation: 

C n e~^ A2 = (VR e-h^ 2 -™** , (i.i) 



where R and A are n x n matrices and C„ is a normalisation factor independent of the matrix 
A. When matrices R and A are, for example, complex Hermitian, the volume element can be 
chosen as T>R cx Yii<j d [R- e ^y] d [Imi?y], and the above integral amounts to a product of 
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standard Gaussian integrals over independent degrees of freedom, the identity ([1.1)1 following 
immediately. The same method works obviously for the real symmetric matrices. On the 
other hand, in these simple cases we also have a freedom to go to "polar" coordinates in the 
standard way. For example, for the complex Hermitian case [10] 

R = U- 1 diag{ Pu ... lPn )U, VR(xdfi H (U)dPA 2 [P], (1.2) 

where U G U(n) is a unitary matrix of eigenvectors, and P = diag(pi, . . . ,p n ) is the real 
diagonal matrix of the associated eigenvalues of R, with dfiH(U) being the corresponding 
invariant Haar measure on the unitary group and A[P] = Oi<j(Pj ~ Pj) standing for the 
Vandermonde determinant factor. Similarly, for the real symmetric matrices 

R = 6- 1 P6, VR(xdfj. H (0)dP\A[P]\, (1.3) 

with O S O(n) being an orthogonal matrix. 

In the problems of interest in electronic transport and random matrix theory the structure 
of the matrices R and A is however restricted by the underlying symmetries of the system, 
and is rather non-trivial, see [11] for a review. For the simplest choice of the disordered 
Hamiltonian corresponding to a system with broken time-reversal symmetry, one of the 
legitimate choices of the integration domain for R is due to Schafer and Wegner[4]: 

R = ATT f + iP, (1.4) 

where the matrices T must be chosen in the pseudounitary group: T G U (711,712). The 
matrices P are Hermitian block-diagonal: P — diag(P ni , P n2 ) — P : an d A > is an arbi- 
trary positive number. For Hamiltonians respecting time-reversal symmetry the integration 
domain R is essentially of the same form, but with matrices P real symmetric block-diagonal 
and the matrices T taken as elements of the real pseudoorthogonal group: T G 0(ni, Ua). 

Although the Schafer- Wegner parameterisation of the integration manifold is correct, an 
accurate verification of the main formula Eq. ljl.ip is not at all trivial, and was provided 
only recently [11]. Actually, this type of parametrization has never been widely used in the 
physical literature. Instead, an alternative parameterisation due to Pruisken and Schafer [5 
has been assumed, tacitly or explicitly, in the vast majority of applications: 

R = f- 1 Pf, VR = d^ H {T)dP l dP 2 A 2 [P]. (1.5) 

Here we assumed the case of broken time-reversal symmetry, T G \J(ni,n 2 ) and P = 
diag(P ni , P n2 ), w hh P ni and P„ 2 being real diagonal, d/xjj(T) being the invariant Haar 
measure on the pseudounitary group and A[P] = Ili<jfe ~ Pj) ^ s the Vandermonde de- 
terminant factor. Apparently, this parametrization is a complete analogue of that in the 
formula (jl.2p . specified for the pseudo-unitary symmetry. 

Similarly, one expects that a natural analogue of (|1.3p for the preserved time-reversal 
Hamiltonians and emerging real-hyperbolic domain should be 

R = f- 1 Pf, VR = d f i H (T)dP 1 dP 2 \A[P]\, (1.6) 

where this time T G 0(rii, n 2 ) is the corresponding pseudo-orthogonal matrices. 

To the best of our knowledge, the validity of the Hubbard-Stratonovich transformation 
with the Pruisken-Schafer choice of the integration domain has not been carefully checked, 
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but rather taken for granted. In fact, the simplest version of the "deformation of contour" 
argument used to verify the transformation for the Schafer-Wegner domain fails for the 
Pruisken-Schafer choice [TTj . and this raised legitimate doubts on its validity in general, see 
also P) . 

Given the widespread use of the Pruisken-Schafer parametrisation, as well as known tech- 
nical advantages of working with it in some microscopic models, the situation clearly calls for 
further analysis. To this end, a rigorous proof of the validity of the Hubbard-Stratonovich 
transformation for the general pseudounitary Pruisken-Schafer domain (|1.5[) was given for the 
first time by one of the authors |13j . In the same paper a variant of the Hubbard-Stratonovich 
transformation for disordered systems with an additional chiral symmetry was also provided. 

On the other hand, the problem of verifying Hubbard-Stratonovich transformation for the 
general real pseudoorthogonal Pruisken-Schafer domain (|1.6[) turned out to be much more 
challenging due to serious technical difficulties to be discussed later on in the text of the 
paper. Only the simplest, yet non-trivial case 0(1, 1) was managed successfully in |13j . and 
we summarize the results of that study below. The integration domain on the right hand 
side of Eq. (jl.ip is given explicitly by 

R = f- 1 Pf, (1.7) 

where 

- / cosh6» sinh6» \ 0(1,1) ~ , 

T= { sinhfl coshfl J 6 0(1) x 0(1) ' ^ P = to*fa>P*)- d" 8 ) 

The matrices A in Eq. (jl.ip has the following form 

A = ( a a 1 ~o 2 ) ' with ai > °' 02 > °' 1°' < V a i°2 ■ (1-9) 

As has been shown in [13] the desirable form (jl.ip of the Hubbard-Stratonovich transforma- 
tion is only possible after one makes the following choice of volume element on the integration 
manifold 

dR = (pi -p 2 )d Pl d P2 d0, (1.10) 
whereas the would-be " natural" choice of the non-negative volume element 

dR = \pi — pi\dp\ dpi dO , 
as in (|1.6p . can not yield a Gaussian function in the left-hand side of (jl.ip . 



In the present paper we continue that study by considering two more specific cases - 
0(2, 1) and 0(2,2), and investigating in detail the validity of the Hubbard-Stratonovich 
transformation for the corresponding real hyperbolic domains. Note that for practical needs 
of the theory of disordered systems 0(2, 2) is the most important case related, in the su- 
persymmetric version, to the basic object of the theory, the so-called two-point correlation 
function of resolvents of the random Schroedinger operator, see e.g. [T| [TT|. 

In both 0(2, 1) and 0(2, 2) cases we are able to show that the naive choice of the measure 
Eq. (|1.6p is never possible, but the Hubbard-Stratonovich transformation (jl.ip can be saved 
provided we make a suitable alternative choice of DP. These examples naturally suggest 
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to put forward the following conjecture on the correct form of the Hubbard-Stratonovich 
transformation on a general 0(m, n)— invariant Pruisken-Schafer domain. Define 

R = f- 1 Pf, P = diag(A,A) = diag(pii,...,pi m) p 2 i,---,P2n) (1.11) 
and the volume element 

m n 

VR = dfi H (T)VP, PP = |A[A]|-|A[P 2 ]| JJ JJ(p u -fty). (1-12) 

i=l i=i 

where |A[P]| is the absolute value of the Vandermonde determinant, and cI[1h(T) stands 
for the invariant measure on 0(m,n). Further assume that the the real matrix A is of the 
form A — A + L, where A + is positive definite and L is the signature matrix L appearing 
in the definition of the pseudoorthogonal group 0(m, nj^. Then the Hubbard-Stratonovich 
transformation over the Pruisken-Schafer type of real hyperbolic domain is given by 

f dii H {f)e-™ f ~ lpfA (1.13) 

•/ O(m.n) 

- lTr-4 2 

= const, e 2 . 

The formula Eq. (|1.13[) is the central message of our work. The crucial difference of 
the choice (|1.12j) from the naive choice of the measure Q1.6P is the absence of modulus for 
the factors rii=i ET^i (Pm ~ P2j)- This forces the volume element to change sign inside 
the integration domain, in contrast to the conventional measures (densities) which are al- 
ways positive as in e.g. Eq. (|1.3p . Such feature does not however in any way invalidate our 
Hubbard-Stratonovich formula, which should be interpreted as follows. The actual sign of 
VR is determined by the inequalities between pi's and P2's. An ordered sequence of the pi's 
and p2's thus defines a sub-domain of R on which the sign of VR is fixed. Without loss of 
generality, we can assume pn > pn > ■ ■ ■ > p\ m and P21 > pn > ■ ■ ■ > pin ■ Then it is 
clear that the domain of integration in R is a union of altogether (m + n)l/m\n\ such disjoint 
sub-domains. Labelling a particular choice of the sub-domain of this sort by D a and defining 
sgn(cr) to be the sign of the volume element VR on D a , the left-hand side of the integration 
formula we discuss is given by 

/ VR e -^-™RA = J- sgn((j) f e -iTrfl 2 - lT rM _ (1 w) 

Interpreting our formula in this way, we always integrate over each sub-domain D a with 
the well-defined positive measures \ VR\, but the l.h.s. of Eq. (|1.13[) is given by an alternating 
sum of integrals on the disjoint sub-domains of R. We believe this coordinated change of 
sign is absolutely necessary to ensure the Gaussian form of the result of the integration, the 
conviction being based on the example of |13] and the results of the current paper. 

We consider verification of this conjecture, as well as the discovery of a general mechanism 
which ensures its validity to be a challenging problem reserved for a future research 0. 



VR 



iTxRA 



VP 



E ! 

i=l 



he pi 

3 = 1 



1 Such matrices can always be brought to a real diagonal form by 0(m, n) rotations, see e.g. Appendix B of 
the paper [15] . 

2 A method of proving the validity of the above conjecture in the general case 0(m, n) has recently been 
proposed by M. R. Zirnbauer and the present authors, and will be published elsewhere |16j. 
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2 Verification of the conjecture for 0(2,1) case 



In this section, we consider the Pruisken-Schafer type of parameterisation of integration do- 
main Eq. Ijl.llj) with T being an element of the real pseudoorthogonal group 0(2, 1). The 
real matrix A in Eq. (|1.13p is assumed to be of the form A = A + L, where A + is positive def- 
inite and L is the signature matrix L = diag(l, 1, —1). As mentioned above, such matrices A 
can be always diagonalised as A — T~ l hT, with T £ 0(2, 1) and A is a real diagonal matrix. 
By exploiting the invariance of the Haar measure we can safely choose A to be diagonal, as 
this choice obviously does not change the result of the integration. 

Implementing the Pruisken-Schafer parametrisation , the integral on the right hand side 
of Eq. (| 1 . 1 3(1 is of the form of 

3 

I%f 1] = [vR e -h^R 2 ^RA = V p e -^R P * f dn{f)e-™ f ~ lpfA , (2.15) 

J J-oo JO(2,l) 

where P = diag(pi,p2,P3) and dp(T) is the invariant Haar measure on 0(2, 1). The crucial 
point is that we have to choose the volume element VP to be, cf. Eq. (|1.12|1 . 

VP = \pi -p 2 |(pi ~P3){P2 - p 3 )dp 1 dp 2 dp 3 . (2-16) 

We are going to demonstrate that it is only this choice that validates the Hubbard-Stratonovich 
transformation for our choice of the hyperbolic domain. 

Note that the integral over the pseudoorthogonal group 0(2, 1) on the right hand side of 
Eq. (|2.15p is of the type of the Harish-Chandra-Itzykson-Zuber integral. Although integrals 
of this type have been known long ago for unitary groups [17] and extended more recently to 
pseudounitary groups [18], their analogues for (pseudo)orthogonal groups, which is relevant 
here, remains largely an open problem in mathematical physics, although a few interesting 
insights were obtained very recently [HI [20] . 



2.1 Particular example of the 0(2,1) Hubbard-Stratonovich trans- 
formation 

To elucidate main points of the calculation we first consider a special choice of the (diagonal) 
matrix A, that is 

A = diag(x, x, z) =► e-? TlA2 = e -^ 2 +*) . (2.17) 

Since AL = diag(x, x, —z) > according to our assumption, we have to require x > > z. 

The calculations will be simpler as such A effectively replaces the integration over the 
whole group 0(2, 1) with one over the non-compact Riemannian symmetric space 

0(2,1) „ S0(2,l) 



0(2) x 0(1) S[0(2) x 0(1)] ' 

Denote dfj,(S) the 0(2, 1) invariant measure on the non-compact Riemannian symmetric 
space G/H, with G = 0(2, 1) and H = 0(2) x 0(1). For our special choice of the matrix A 
we obviously have 

dp(f ) e-^rf-ipfA = I M§) , <t- f■<^ ( (2 19) 

0(2,1) JG/H 
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so that Eq. (|2.15j) assumes the following form 

3 

PPe^-/ d^S)e- iTrS ~ lpSA . (2.20) 

-oo JG/H 



To perform the integration over the coset space G/H it is convenient to parametrise G/H 
with the projective coordinates (Z, Z T ). To this end, we introduce a 2 x 1 real matrix Z as 

Z =( Zl ] with the constraint 1 - Z T Z > , (2.21) 



Z2 

in terms of which the matrices S on G/H are given by 



a_i {i-zz T y* z(i-z T zy^ . 
b 1 ^(l-^)-l (i-z r z)-* ' • 



It is direct to check that S~ 1 {Z, Z T ) = S(-Z,-Z T ). The invariant measure d/i(S) in 
projective coordinates can be calculated in the standard way[21] and is given by 

, dZdZ T . 

MS) = JT^W (2 ' 23) 

where dZdZ T = dz\dz 2 and the integration domain is as specified in (|2.21|) . Make the 
following change of variables 

I Zl=rCOS0 a , r G [0,1] and0G [0,2tt]. (2.24) 
I z% = r sin 

The integration on the right hand side of Eq. (|2.19p can be written as 

f 1 rdr f 2v i f r 2 x — Z 

/ T, oTT / d6,ex Po1-j 2^ ~ ~^ 2 ) cos2d + i a^i+P 2-2 ^) 

io (l-r)2 Jo ^L 1 " r i — r 

+ [a?(pi + jj 2 + 2 P3 ) + z(pi + pa)] 1 • (2.25) 



The integral over 8 yields the standard Bessel functions in view of J* d<p e ll3cos ^ — ttJq(P) , 
and introducing a new variable t = * » , we rewrite (|2 . 25[) as 



e i(x( Pl +p 2 )+zp 3 ) / J 

Jo VTTt 



-(x- z){px -p 2 ) 



e ^(x-z)(p l +p2-2p 3 ) ^ (2.26) 



Now we need to substitute Eq. ()2.26|) into the right hand side of Eq. (|2.20j) and to integrate 
over P, that is 

= / -== I TIP expj - J > >J + + pa) + SP3) 



r7m/_> pexp {-5g , '' +i(i& " 



it 



-(x-z)(pi -p 2 ) 



(2.27) 
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After a straightforward, but lengthy calculation we arrive at the following result 

I H r } = ^Fl(*-z) 2 ]e-^ 2+ * 2 \ (2.28) 

where 

f°° dt / 1 \ 

F(a)= cxp ( --ft 2 +t)a ) [l - a(2t 2 + 3t + 1)] . (2.29) 

Note that the expression Eq. (|2.28[) contains already the Gaussian factor of precisely the 
form required by (|2 . 1T|) . Unfortunately, that factor is multiplied with a function F[{x — 
z) 2 } dependent on the combination a = (x — z) 2 , the fact seemingly incompatible with the 
Hubbard-Stratonovich transformation. Miraculously enough, this factor is an a— independent 
constant! To verify this, we define y = y/1 + 1, and carry out the integral explicitly: 

r°° dt ( 1 \ 

F(a) = J -j== cxp (--(t 2 + t)aj [l - a(2t 2 + 3t + 1)] 



/OO 
dy e^(-^~y 2 )) [l - a(2y 4 - y 2 )] 



= 1 - lim ycxp(-^^ % = 1 . (2.30) 

At the last step, we used the fact that a is strictly positive, as the case a = should be 
excluded from the very beginning. Indeed, a = implies x — z, contradicting to the original 
requirement x > > z. 

2.2 General calculation for 0(2,1) case 

Now we are ready to present the complete proof of the Hubbard-Stratonovich transformation 
over 0(2,1) domain. In the general case we have A = diag(a;i, X2, z) = A\ + A2 where 
Ai = diag(x, x, z) is the part considered in the previous example, and A2 — diag(w, — w, 0). 
Here we defined the variables x = (xi + X2)/2, w — (x% — X2)/2. Our starting point is again 
Eq. (|2.15p . but we now have 

I°f 1] = IvR e -^-i^A 



= [°° VPe i S 1 P ' ( e- iTrS ~ lp§Al dn(S) [ d^H) e- iTrS ~ lp ^ kA ^\ 
J -00 Jg/h Jh 

(2.31) 

where we assume G = 0(2, 1), H = 0(2) x 0(1) and S = G/H as before. 

The integration over H goes effectively over the group SO (2) and the corresponding 

matrices can be parametrized in a standard way as H = ( C0S , S ^ n , | . Using the 
y J y — sin (/) cos tp J 6 

same parameters for the coset matrices S as in the previous section, we then find 

Tr S~ 1 PshA 2 H- 1 = A cos 20 + B sin 20 , (2.32) 
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where 



W _. i [(1 + v 7 ! - r 2 ) 2 + 2r 2 cos 20 + cos46»(l - v 7 ! - r 2 ) 2 ! 
1 — H) L 



Pi 



+ [2r 2 cos26 1 - (1 + Vl-r 2 ) 2 - cos 40(1 - \f\ - r 2 ) 2 ]p 2 - 4r 2 cos26»p 3 

B = —. ~ W „. | [2r 2 sin 26> + sin 46>(1 - \/l - r 2 ) 2 lpi + [2r 2 sin 29 - sin 40(1 - v 7 !" 
4(1 — r z ) y 

- Ar 2 cos 28p 3 



The integration over <fi is easily performed according to the formula 
so that 



>exp U cos ( 



/ .///■://:, ' n< > " = J (v^+~B2") 

is 



r 2 ) 2 ]p 2 
(2.33) 

(2.34) 
(2.35) 



This should be inserted into Eq. (|2.31j) . and remembering Eq. f|2.25j) - ()2.27[> . we arrive at 



I 



0(2,1) 
HS 



P OO 7^ />OC />27T rf' X 

Jo J VP Jo d9 I 2 + + K) + ZP3) 



+ |(x - z)( Pl +p 2 - 2p 3 ) + |(x - z)(pi - pa) cos 6>j J ( V^ 2 + £ 2 ), 

(2.36) 

where again XJ>P is given by Eq 1|2.16|) . 

Note that variable 'w' responsible for the difference from the example considered in the 
previous section enters the formula only via the combination \J A 2 + B 2 . A way of evaluating 
the above integral for w ^ is to expand the Bessel function in Taylor series with the n-th 
term proportional to w 2n , to integrate each term separately, and then re-sum the series. A 
straightforward implementation of this program is however not immediate, and necessary 
steps of the proof are given in ApplAl where it is shown that 



it 



rO(2,l) 
L HS 



const exp 



x — w — 



const exp 



/ 2 i 2 i 2\ 

-( Xl +x 2 +z ) 



(2.37) 



in precise agreement with the structure required by the Hubbard-Stratonovich transforma- 
tion. 



To summarize, we have demonstrated that for any A = Tq diag(a;i, z) T 1 and To G 
0(2, 1) holds the identity 

J VR e -^R 2 -^RA = congt e -lTri 2 t (2 3g) 

provided the volume element VP for the P integral is chosen in accordance with Eq. (|2.16|) . 



For the sake of comparison, one may try to repeat the above calculation with the "naive" 
choice of measure DP = |A(P)|n- = i^ instead of Eq. (j2~16j) . We show in App. |B] that 
such a choice invalidates the Hubbard-Stratonovich transformation. As another compar- 
ison, we also provide similar calculations in App. [C] for the compact counterpart of this 
Pruisken-Schaier domain corresponding to the group 0(3). 
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3 Results for the 0(2,2) case 



In this section, we carry out the detailed calculation for the Hubbard-Stratonovich trans- 
formation over the 0(2,2) Pruisken-Schafer domain. As the calculation turns out to be 
quite technically cumbersome, we restrict ourselves with the simplest non-trivial choice 
A = diag(x, x, z, z), with x > > z. Consequently, the integration domain T — 0(2,2) 
effectively reduces to the non-compact Riemannian symmetric space (coset space) 

0(2,2) ^ SO(2,2) 
0(2) x 0(2) S[0(2) x 0(2)] ' [ ' 

Parameterisation of G/H, where G = S0(2, 2) and H = S[0(2) x 0(2)], with the projective 
coordinates Z and Z T is again in the form of Eq. (|2.22|) with Z and Z T being real 2x2 
matrices chosen in a way ensuring that the matrix 1 — Z T Z is positive definite: 



z\ z-i 

Z 3 Z4 



with 1 - Z T Z > . 



(3.40) 



We aim to prove the validity of the Hubbard-Stratonovich transformation with the 
Pruisken-Schafer parameterisation Eq. (|1.5|) . where Te 0(2, 2) and P = diag(pi,p2,P3,P4)- 
To this end, we need to demonstrate that the following integral 



0(2,2) 
HS 



00 . 4EpI 

VP e 

— oo 

VPe a £ P ' 



d(j,(T)e 



-iTvT~ PTA 



0(2,2) 



G/H 



d»(S) e - m§ ~ lpSA (3.41) 



is, up to a constant factor, a product of Gaussian factors. The invariant measure dfi(S) here 
is calculated in the standard way and is equal to |21j 

dZdZ T 



dfi(S) = 



det(l - Z T Zf 



(3.42) 



where dZdZ T = dz-ydzidz^dz^ 



To carry out the integration over the coset space we introduce the polar coordinates 
parametrization for real matrices Z. This amounts to diagonalizing Z by two orthogonal 
rotations as 



Z = O x 



r 
s 



2 , where r, s G (-oo, oo), Oi,0 2 eSO(2). 



(3.43) 



A standard calculation (App. |D]) shows that the Jacobian induced by changing variables 
from Z, Z T to the polar coordinates is simply |r 2 — s 2 \. We have accordingly 



dZdZ 1 



\r 2 - s 2 \drdsdn{0 1 )dn(0 2 ), 



(3.44) 



where dpk(0\) and d/i(Oa) are the invariant Haar measure of SO (2). Using the polar coordi- 
nates the integral over coset space takes the form 



G/H 

D(r,s) 



d^{0\) exp < iTr 



SO(2) 



d/i(02) exp < iTi 



SO(2) 




(3.45) 
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where we denoted D(r,s) = \r l — s 2 \drds/(l — r 2 ) 2 (l — s 2 ) 2 



The two integrals over 0(2) group manifold in Eq. (|3.45j) are easily carried out using the 
formula 



/ d^(0)exp liTr O 



01 \ J / b X 

a 2 J b 2 



exp 



-(ai + a 2 )(bi + b 2 ) 



Jo 



(ai - a 2 )(bi - b 2 ) 



(3.46) 



Introducing at the next step the variables u = * 2 and v = 5 , we rewrite the resulting 
integral in Eq. (|3 .45[) as 



(p 3 +p 4 ) +iz (p i +p 2 ) 



\u — v\dudv ( i 

r=, \ \ ex P \ t;( x ^ z ){Pi+P2-P3-P4)(u + v) 

1 i/u(u - l)yjv(v - 1) I 2 



Jo 



-(x- z)(pi -p 2 )(u- u) 



Jo 



-(X - z)(jJ 3 -P4)(W - I') 



(3.47) 



Now we have to perform the integration over variables in P. As in the previous section, the 
crucial point is to choose the volume element DP in accordance with our main conjecture, 
that is 

4 

DP = \pi -P2KP1 ~P3)(pi -Pi)(P2 -Pz){P2 ~Pa)\pz -Pi\\\dp l . (3.48) 

i=l 

The remaining steps are lengthy but straightforward. After a few variable changes we arrive 
at 



,-SO(2,2) _ _k_ 

HS ~ 128" 



= T^:F[a] exp [- 



2 21 

X — Z \ 



with a = x — z and the function T[a] given in terms of a double integral as 



(3.49) 



dt e 



_^ r(t-i) 2 \aH\t 2 - v) - aH 2 + 1 



y/[(t +1) 2 -v][(t- I) 2 -v] 
Integrating over v and defining a new variable x — ^4p, we get 



e 4 dv 



(3.50) 



7T 

128 J 1 



T[a} = — I dxe-^ 2 -^^{[a 2 (2x-l) 2 -2] 2 Z 1 [l,t^(—y,-a 2 (x-l) 2 ] 



1 3 ,x — l x 2 
\a\x - lf(2x - l) 2 $i [2, \, | (^-l) 2 , -« 2 (x - I) 2 ] 



(3.51) 



where $i is the degenerate hypergeometric series of two variables defined as 



$i[a,/3,7,x,y] = ^2 



m rn n n 
X y 



^ n (7)m+n ™!n! 
m,n— v 7 



(3.52) 



From Eq. (|3.49|) we see that only if the factor J- [a] is independent of its argument a = x—z 
the whole expression I^fg 2 ' 2 ^ can be in the desired Gaussian form. It needs only a few lines 
of Maple or Mathematica code to check numerically that actually F[a\ = 1, see FigfTJ 
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F(a) 



1 

0.8 
0.6 
0.4 
0.2 

20 40 60 80 100 a 

Figure 1: Function F(a) = 1 is a constant which does not depend on a. 



Unfortunately, we were not able to find a way of verifying this miraculous identity ana- 
lytically, as we managed to do in the previous case of 0(2, 1) integral. Nevertheless, we do 
not think the numerical data leave any doubt in the validity of our claim. 

In conclusion, the above calculation shows that for A = To diag(ir, x, z, z) Tq 1 and Tq £ 
0(2,2), 

j VR e-^« 2 -™ = const e -^ A \ (3.53) 

provided measure for P integral is chosen to be Eq. (|3.48|) . 

It is again interesting to check what will be the result if we choose dP = |A(F)| J]Li d Pi 
instead of Eq. (|3.48p . It is shown in App.lElthat this choice will make the Hubbard-Stratonovich 
identity invalid. 
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A Proof of Eq. ( IQ7I) 

Introducing the set of new variables 

a=\(pi-P2), b=^{ Pl +p 2 ), c = p 3 (A.54) 
and defining t — jzr~s, we can rewrite Eq. (|2 . 36[) as 

r°° dt r°° r°° ( r 2 

r O(2,l) _ / / ^ / ^ _ c)2 _ fl 2j cxp I _ a 2 _ b 2 _ _ + ^ x{b + c j + 

J vl+tJo J-oo I 2 

+ i cos26» ta(x -z) + i(l + t)(x - z)(b - c) \j a (^A 2 + B 2 ) , (A.55) 



HS 
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w 



where we have used Eq. (|2. 16|) . Next, we verify that 

A 2 + B 2 = w 2 {t 2 (b- c) 2 + 2t(t + 2) cos 29 a(b - c) + a 2 [t 2 cos 2 29 4( / + 1 i| 

[t(b -c) + a(t + 2) cos 26>] 2 + a 2 [4(t + 1) sin 2 29] 
= C 2 + D 2 , (A.56) 
where we defined 

C =w[t(b - c) + a(t + 2) cos 29] 



D =2wa^/t + 1 sin 29 . (A.57) 
This allows us to write (cf. fl2~34]) ) 

Jolv 7 ^ 2 + S 2 ) = - / d0exp(icos0C + isin0D) . (A.58) 



Using this representation of the Bessel function in Eq. (|A.55p and defining y = 1/ vl + f , we 
can readily carry out integrals over a, 6, c and 0, and get 



= exp[-x 2 - y ]F(«;) , (A.59) 



r O(2,l) _ 

where we defined 

- [iu 2 - (a; - z) 2 ]y 2 - [w cos <?(>(y 2 - 1) + (a; - z)y 2 ] 2 \ . (A.60) 



/>00 />7T 

F(w) — dy d(f>exp 
J i Jo 



Note that although the above integral formally seems to depend on both w and (x — z), 
we shall see below that it is a function of w only and is actually independent of the second 
combination. 

To calculate F(w) we find it convenient to apply first the standard Hubbard-Stratonovich 
transformation and " linearise" the second term in the exponent by introducing an auxiliary 
Gaussian integral: 



oo />7r />oc 



F(w) = I dy J d(/) I d/iexp<{ ~ [w 2 - (x - z) 2 ]y 2 - h 2 



+ 2ih [w cos cf)(y 2 - 1) + (x - z)y 2 ] j . (A.61) 

Integration over (f> yields the Bessel function which can be expanded in its Taylor series, and 
the Gaussian integral over h can be performed. In this way we find 



00 roc f 

F(w) = const J2 w2n / d y e- {x - z)2{yi - y2) i 



\-2{l-2y 2 ){x-zf 



C n + 2(1 - 2y 2 )C n - 1 
(A.62) 

where we defined for n > 

in — my. mlm\ 4 k k! 2m — 2k )! 

m=0 v ; fc=0 v ; 
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and C n = for n < 0. In particular, the definition above implies 

_ (-)" 



(A.64) 

y=l 



F(w) can be found as we are now able to perform the integrations on the right-hand side of 
Eq. (jX62|) as 

{i 
[- - 2(1 - 2y 2 )(x - z) 2 ]C n + 2(1 - 2y 2 )C„-i 
y 

n-l 

[C n + y 2 (^l)^fl, l! ,C 1 ] 



i I V 

e 



-(x-z)\y*-y 2 ) 



(-)" 



71! 

Here a n ^s are coefficients satisfying the following recursive relations 



(A.65) 



2 1 

a n , n -i = -, and a n , t = a n -x,i, i = 0,...,n— 2, a 10 = 2. (A. 66) 

n n 

In the last step of Eq. (|A.65[) we used the fact x — z > 0. We finally see that Eq. (|A.65|) 
implies the desired Gaussian expression 

F(w) = const e~ w2 (A.67) 

Finally, substituting F{w) oc e~ w2 back to Eq. (|A.59|) completes our proof of Eq. (|2.37[) . 

B Calculation with the naive choice of the volume ele- 
ment DP for 0(2, 1) case 

In this appendix, we show the Hubbard-Stratonovich transformation for the 0(2, 1) Pruisken-Schafer 
domain is invalid if the volume element is chosen to be DP = |A(P)| Yli=i ^Pi- 
Starting from Eq. (|2.27p . we make a change of integration variables as in Eq. (|A.54|) . 
Then we write Eq. ([2~27]) as 



Z°^=/ -^=\ ppe-^-^-^WMMj^t^.^], (b.68) 



x dt_ 
o VT+tJ- c 

where 

DP = 2\a((b-c) 2 -a 2 )\dadbdc. (B.69) 

We rewrite the above integral as 

t O(2,i) _ t O(2,i) r O(2,i) m 7rn 

where we defined 

rOC 7j. poo 2 

jOg.i) = y 2| a |((6-c) 2 -a 2 )dad6dce- a2 - b2 -^+ l ( 2a;b + 2C )+ Jt ( a; -^( fc - c ) 

Jo [*(a: - z)a] (B.71) 
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and 



r O(2,l) 
'HS,2 



dt 



o y/T+tJo 



b-c 



4\a\(a 2 - (b-c) 2 ) da 
Jq [t(x — z)a] . 



dbdc 



-a 2 -b 2 -^ T +i(2xb+zc)+it(x-z)(b-c) 

(B.72) 



Let us stress that it is the contribution 2 which encapsulates the difference between 
the definition DP = |A(P)| n»=i dp% which is positive definite and Eq. (12.16|) which is sign 
indefinite. Such a term has cancelled out when we the volume element was chosen to be 
Eq. (|2~T6f . The first contribut ion I H g ^ ^ is nothing else but the I H g ' calculated in the 
Section 2, and we proved it is in the Gaussian form 



rO(2,l) 
*HS 



const, exp 



-(2x 2 + z 2 ) 



r O(2,l) 
-HS.2 



(B.73) 



In the remaining part of this appendix we will demonstrate that the first contribution 
2^g2 i s n °t m a Gaussian form, thus invalidating the Hubbard- Stratonovich transformation. 



Define m — b + c, n = b — c and integrate over m. We get 
0(2,1) _ 2\/67r 

dt 



-HS.2 



3 

/ da 

Jo jo 

r O(2,l 



1 Z" 00 i f 1 

exp [ — -(2x + z) 2 ] / dn n A exp < — (- 



- P + a 2 )n 2 + in(x - z){- + t) 



V1 + * 



a(a — 1) Jq 



t(x — z)na 



(B.74) 



It is clear that 1 HS 2 will be in the desired Gaussian form if integral part of the above 
formula is oc exp(— (x — z) 2 /3). To check this, it is sufficient to consider a special case x — > z, 
i.e. \x — z\ -C 1. In this limit, we can approximate the integral by setting the argument of 
the Bessel function in the integrand to zero. This gives 



r O(2,i) 

'HS,2 



The integral over a is simply 



dt 



y/T+t 



1 2 

dn n A exp < n 2 + in(x — z)( — h t) 



daa(a — 1) exp(— n a ) 



(B.75) 



f 1 1 

/ daa(l- a 2 )exp(-n 2 a 2 ) = — ? [exp(-n 2 ) + n 2 - 1]. (B.76) 



Carrying out the standard Gaussian integrals over n, we get 



t O(2,1) -< 2 " + '' 



Vl + i 



- [(z - z) 2 (3t + 2) 2 - 2]e-&<*-' )2(3t+2)2 



e -3s(*-*) 2 (3t+2) 2 



(B.77) 



The integral over t is divergent if x — z = 0, as expected, and in the limit \x — z\ <C 1 it is a 
well-defined expression dominated by t ~ (a; — ^ 1 so that ^go*' ~ C 35 — z) -1 / 2 . Such 
a pre-exponential factor clearly precludes the expression to be in the desired Gaussian form. 
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C Calculations for the standard case 0(3) 



In this appendix, we repeat calculations similar to those in section 2 and App. |5J but this 
time for the compact case of 0(3) group. Although the Hubbard-Stratonovich transforma- 
tion for 0(3) symmetry is trivially valid in the original formulation, it is instructive to have 
a comparison between 0(3) and 0(2, 1) in the polar representation, as it helps to understand 
peculiarities of the non-compact case. 

First, we consider an integral similar to Eq. (|2.15|) . with integration of T going this time 
over 0(3) instead of 0(2,1). We consider only the simplest case setting A — diag(ir, x, z) 
and deal with the following integral 



3 

I°f = [vR e -i™ 2 -™ M = f°° VPe k ^ A [ d»(S) e^-PSi ; (CJ8) 

J J-oo Jg/H 



Ig/h 

where G = 0(3) and H = 0(2) x 0(1). Elements of this compact coset is parametrised as 



_/ (i + zz T yi z(i + z T zyi > 

S -9*-\ Z T(l + ZZT)-i (l + ZT Z )-i >' (CJ9) 



where we introduced the 2x1 real matrix Z as the convenient coordinate on G/H , with 

Z = ^ Zl ^ , with zi and z 2 arbitrary real. (C.80) 

Similar to the non-compact case, s _1 (Z, Z T ) = s(—Z, —Z T ). The invariant measure dfi(S) 
in the projective coordinates is given by 

dn{S) = — F , where dZdZ T = dz x dz 2 . (C.81) 

(1 -f- Z^ Z^ 2 

The integration over the coset is now straightforward and calculations are done parallel to 
those in section 2. After some algebra and a few changes of variables, we get 



.0(3) 
1 HS 



f 1 dt f°° ( 1 3 

Jo V^t J VP 6XP I 2 P * + i{X ( Pl + P2) + ZP3) 



it 

+ -^( z ~ x )(Pl +P2 ~ 2p 3 ) } J 



-(x-z)(px -p 2 ) 



(C.82) 



The difference between Eq. (jC.82l) and Eq. (|2 . 27[) is due to the difference between compact 
and non-compact integration manifolds. 

A crucial difference in the 0(3) case is that the volume elements VP in the above formula 
is VP = |A(P)| J]Li d Pi> instead of Eq. (|2~T6| . We have seen in App. Ethat this choice of 
VP when applied for 0(2,1) symmetry would yield a form which is not Gaussian. In the 
remaining part of this appendix we show that in the case of 0(3) the result is in contrast 
Gaussian. 

Define the same set of integration variables as Eq. (IA.54P and use them in Eq. (IC.82|) . 
We have 

T o(3) _ T o(3) T o(3) (r 8 „, 
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where we defined 

J°f 1= / -fL= / 2\a\{{b~ c) 2 - a 2 )dadbdc e' a2 - b2 -^ 2xa+ ^- lt ^^ a -^ 
Jo y/l — t J-oo 

Jo [t(x - z)a] (C.84) 

and 

,1 ^ f \b-c\ 
'o yl—t Jo 

J [t(x - z)a] . (C.85) 

Note again that T^g \ corresponds to the definition Eq. (|2.16[) and Z^S 2 emer g es only be- 
cause the volume element is positive definite in the current case. 



jOW = J J 4 | a | (a 2 _ (6 _ c)2) da / ^ g-.'-*- V44(2xa+*c)-tt(*-,)(«-c) 



find 



First, we deal with Carrying out simple Gaussian integrations over a, b and c we 



« = ^i(^- (2 * 2+ * 2) , (C86) 

where 

T x {x, z) = J -^L= exp (^-i(f 2 + t)(x - zf^j [l-(x- z) 2 (2t 2 + 3t + 1)] . (C.87) 

Using a — (x — z) 2 and y — y/1 — t we immediately see that 

F x {a) = f l -jt= exp (-\{t 2 + t)(x z) 2 ) [l-(x- z) 2 {2t 2 -3t + 1)] 

dy exp(-^(y 4 - y 2 )) [l - a(2y i - y 2 )} 

= 1 - lim y exp(- ayl( y 2 Z jj ) = i . (C.88) 

Here, the integral over y is the same as the one in Eq. (|2.30[) with different upper and lower 
limits, but the result is the same. This completes our proof that \ indeed in the Gaussian 
form. We also note there is certain kind of duality between I^s l an0 - ^hs ■ 

As we know already, the integral l& °f the Gaussian form cx exp[— x 2 — z 2 /2], and 

we have just shown that the same holds for Tf^\, the second terms X H ^\ can then only 

be either or the same Gaussian form as I^g \ ■ To see which is the case, it is sufficient to 
consider the same limit x — > z as we did in App. [B] In the limit \x — z\ <C 1, we find 

X^ 3) 2 =/^=/ A\a\(a 2 -(b-c) 2 )da / dbdc e -a 2 -b 2 -^+«2 Xa+zc ) (C 89) 

Jo V 1 — t Jo J-oo 

One can perform all the integrations in this formula explicitly, and show that 

J°g ) 2 cxexp{-i(2a; 2 + z 2 )} (C.90) 

as expected. 
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D Jacobian of the transformation from Z to polar coor- 
dinates 



Write the polar coordinates decomposition in Eq. (I3.43j) as Z = 0\A0 2l where A = 
We have 

dZ = dOiA0 2 + 1 dA0 2 + OiAd0 2 
dZ T = dOlAO Y + O^dAOx + 2 r Ad0 1 ' 

Following the standard way of derivation, see e.g. [21], we have 

d 2 S = Tr dZdZ T = Tr jof d0 1 A0 2 dOl A + AdA0 2 dOl + A 2 d0 2 d6l 
+ OldOiAdA + d 2 A + AdAd0 2 Ol 
+ A 2 dOld0 1 + AdAdO^Ox + Ad0 2 O^AdOjO! 



r 
s 



(D.91) 



Next we define 



50i = 0\dO x 



S0 2 = d0 2 Oi ' 

Recalling that SO\ and 50 2 are skew-symmetric matrices, which can be written as 



SO hl2 
-50 hl2 



so 2 



S0 2 , u 
-50 2 ,i2 



We find 



d z S = Tr<j d 2 A - A^SOiSOx - A 2 S0 2 S0 2 - 2AS0 1 A50 2 

- 2 2rs \ 



( r 2 + s* 



{50i. i 2l S0 2i i 2 ,dr, ds) 



1 



V 



1 / 



/ so 1A2 \ 
so 2A2 

dr 

V da J 



(D.92) 



(D.93) 



(D.94) 



= dx^ijdx 3 . (D.95) 
In the last step the summation over repeated indices is assumed. Jacobian is then given by 

Jacobian = y det g = \r 2 — s \ . (D.96) 



E Calculation with the alternative volume element DP 

for 0(2, 2) case 

In this appendix, we calculate Eq. (13.4111 with the volume element DP = |A[P]| Yit=i dpi 
used instead of Eq. (|3.48[) . We show that by this choice the final result is not in the Gaussian 
form, hence the corresponding Hubbard-Stratonovich transformation can not be valid. 



First we redefine the integration variables 

: Upi +P2) 
: 2 {Pi ~Pi) 
: 2(^3 +Pi) 
■ 2 {P3 - Pi) 



(E.97) 
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Then we have 

|A[P]| = A\bd\-\[{a-c + df -b 2 ][(a-c-d) 2 - b 2 }\. (E.98) 
Use Eq. (|3.47p and the DP defined above to find 

jO(2,2) = f°° \U-V\dudv f D p e -a 2 -b 2 -c 2 -d 2 +i(x-z)(a-c)(u+v) 

HS Jl l)y/v(v— 1) J 

J [b(x -z)(u- v)] J [d(sc -z)(u- v)} . (E.99) 
As in App. [B]we can split into two parts, 

t O(2,2) _ r O(2,2) t O(2,2) ,„ , „ m 

Here, the contribution iSg'i ' s precisely the Ijjg we calculated in section 3, and we 
know it is in a Gaussian form. It is the second contribution, 2^S2^' which arises from the 
difference between the two definitions of DP, and it is given by 



-oo 
oo 



l— c\4-d 



X H [ S? « / dadc I dd I db bd[(a -c + d) 2 - b 2 } [(a - c - df - 6 21 



||o-c|-d| 
u — v I dudv 



^u(u - l)\Jv{v - 1) 
J [b(x- z)(u-v)]J [d{x- z)(u-v)} . (E.101) 

In the remaining part of this section we demonstrate that T^g'^ is not in the Gaussian 

0(2 2) 

form, thus I HS ' is not either. Which means that Hubbard-Stratonovich transformation 
fails with this different choice of DP. 

First, we define m — a + c and n = a — c. It is clear that the integral over m is decoupled 
from other integrations and can be easily performed. Again, it is sufficient to consider the 
limit |a; — ^| <C 1. For the same reason as in App. [Bj we set the two Bessel terms to be 1 in 
this limit. We then have 



I poo poo rn+d 

X H [ sf <xexp{--(a; + z) 2 } / dn j dd j db bd[(n + df - b 2 ][(n - d) 2 - b 2 ] 



n—d\ 

| It — v\dudv f 2 j,2 



exp {-n 2 - b 2 - d 2 } cos [n(x — z)(u + v — 1)] 



yju{u - l)y/v(v - 1) 

(E.102) 

The integral part of the above formula needs to be oc exp{— \{x— z) 2 } in order to makeX^ 2 ^ 
be Gaussian. The remaining calculations are lengthy but direct. We perform Gaussian type 
integrals over b, d and n then define new integration variables X = u + v — 1 and Y = u — v. 
After integrating over Y we get 

/' OC 1 2 2 n V 

X HS2 )(X dX\n(2X + l){8~2a 2 X 2 + y /^ e - £ -^aX{a 2 X 2 -6)Erfi[— ]} , (E.103) 

where we defined a = x — z, and Erfi stands for the error function of imaginary argument. 
Integrating by parts we bring the above integral to the form 

f^2XTl(^-^-^(a^-2)E r fi[f ]). (E.104) 
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Again, in the limit \x—z\ -C 1 the integral over X is dominated by the region X ~ (x— z) 1 3> 
1. Changing variable aX — > X and expanding in terms of \x — z\, we find to the lowest order, 
^-HS2^ = c o — N ~ z l c i + ^((^ — z ) 2 ); with Co and c\ being some constants. In this way one 
finds that * s a function of (x — z), but is clearly not in the Gaussian form. So we 

conclude that T^g is not Gaussian. 
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